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Abstract 

We derive solutions to the Schwinger-Dyson equations on the Closed-Time-Path 
for a scalar field in the limit where backreaction is neglected. In Wigner space, 
the two-point Wightman functions have the curious property that the equilibrium 
component has a finite width, while the out-of equilibrium component has zero 
width. This feature is confirmed in a numerical simulation for scalar field theory 
with quartic interactions. When substituting these solutions into the collision 
term, we observe that an expansion including terms of all orders in gradients leads 
to an effective finite-width. Besides, we observe no breakdown of perturbation 
theory, that is sometimes associated with pinch singularities. The effective width 
is identical with the width of the equilibrium component. Therefore, reconciliation 
between the zero-width behaviour and the usual notion in kinetic theory, that 
the out-of-equilibrium contributions have a finite width as well, is achieved. This 
result may also be viewed as a generalisation of the fiuctuation-dissipation relation 
to out-of-equilibrium systems with negligible backreaction. 



1 Introduction 



Observables in finite-density systems are typically derived from expectation values of 
operators that are not time-ordered. In equilibrium field theory, these expectation val- 
ues can be obtained from analytic continuation of the corresponding expectation values 
calculated in the Matsubara imaginary time formalism. Alternatively, one may calculate 
the expectation values directly in real time. For this purpose, diagrammatic rules can be 
formulated that are known as the Closed-Time-Path (CTP) formalism [H, |2[. For out-of- 
equilibrium systems, the Matsubara formalism does not apply and the CTP formalism 
is the standard method for a description derived from first principles of Quantum Field 
Theory 

While a canonical formulation is also possible jsj, the CTP Feynman rules can easily 
be derived from a CTP generating functional jsf. Moreover, a particularly useful variant 
of the functional approach is to use the two-particle-irreducible (2PI) effective action in 
order to derive Schwinger-Dyson equations. These are self-consistent equations for the 
dressed propagators A(a;, y) and in principle encompass the full real time evolution of 
the quantum system up to two-loop order in the effective action 0, [sj . 

Analytical approaches typically base on linear approximations to the Schwinger- 
Dyson equations. The linearisation may be realised by neglecting the functional depen- 
dence of the self-energy Il{x, y) on the propagators A(a;, y), which implies the neglection 
of backreaction. Moreover, it is often assumed that Il{x, y) is invariant under space- 
time translations, as it is appropriate for a homogeneous, time-independent background. 
The linearisation can therefore be justified by the gradient expansion 0, M, H or by 
the notion of a bath, that does not mediate a backreaction to the out-of-equilibrium 



test field |llHl4j. These approximations are often applicable to weakly coupled systems 
close-to-equilibrium js], [l5| . 

One could therefore argue that instead of the CTP-formalism, one might as well use a 
heuristic and more intuitive Boltzmann approach in order to formulate kinetic equations 
that can be solved analytically, or at least semi-analytically. Moreover, such Boltzmann 
equations can straightforwardly be modified in order to account for a finite width of 
the interacting quasi-particles and in order to account for Fermi-Dirac or Bose-Einstein 
statistics. For finite density QCD, a powerful kinetic theory that includes off-shell effects 



is formulated this way in Refs. |16H21 



The present work is motivated by a number of recent a 



12 



lications of the CTP for- 
In single-flavour 



22-31 



malism to leptogenesis and electroweak baryogenesis 
leptogenesis, the CTP formulation naturally removes the non-unitary overcounting of 
real intermediate states, that has to be performed explicitly when using the heuristic 
Boltzmann ansatz. Moreover, the CTP-formalism readily provides a framework for the 



description of flavour coherence. It can be successfully applied to flavour leptogenesis [28 



and to models of electroweak baryogenesis from flavour mixing (22|, |23|, |3lJ . The approach 



bears the promise of further progress in the directions of A^2-leptogenesis, leptogenesis in 
the weak washout regime and resonant leptogenesis. 

These applications show clear advantages over the conventional Boltzmann approach 
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and justify a further development of the kinetic theory based on the CTP formahsm. 
The improvements from the apphcation of CTP methods to leptogenesis beyond the 
Boltzmann approach additionally encompass thermal corrections in the cut propagators 
of intermediate on-shell particles 12|, |2J-|27| and the treatment of flavour coherence 



28 



It is then natural to calculate corrections beyond these leading order approximations. 
These can originate from higher order derivatives 29|] or from considering the effect of 



the finite width of the singlet neutrino, that decays out of equilibrium [12|, [14 



In Refs. 11-14 



the Schwinger-Dyson equations in the approximation of negligible 
backreaction are solved in order to obtain the propagator of an out-of-equilibrium field 
that includes finite-width corrections. We emphasise that the assumption of a stationary 
bath is equivalent to any approximation scheme where the dependence of the self-energy 
of the out-of-equilibrium field on the average time and on derivatives with respect to 
average time are neglected, such as some variants of the gradient expansion. In the 
case of leptogenesis, the out-of-equilibrium field is a heavy singlet neutrino and the 
background bath consists of the O(IOO) Standard Model degrees of freedom that are in 
equilibrium during the time of leptogenesis. 

Neglecting the backreaction leads to a somewhat curious form of the Wightman 
functions for the out-of equilibrium field [cf. Eqs. fll8|2ip below]. (The Wightman 
functions are the two-point correlations without time ordering. They contain statistical 
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First, there is a 



information about the system, cf. Appendix |S] and Refs. 
contribution corresponding to a quasi-particle distribution that is in equilibrium with 
the background bath. This distribution has a finite width around the quasi-particle pole, 
which is given by the relaxation rate of the system toward equilibrium, in accordance 
with the fluctuation-dissipation theorem for equilibrium systems. In coordinate space, 
the finite width can be identified with an exponential decay in the relative coordinate of 
the Wightman function iA^'^{x,y). Second, there is an out-of equilibrium contribution 
that exhibits no decay in the relative coordinate, or equivalently zero width. 

There is an interesting connection of the solutions [Eq. (fT8|) below], that decompose 
into a finite-width equilibrium and a zero width out-of equilibrium contribution, with 
the problem of pinch singularities in out-of-equilibrium field theory in the real time 
formalism. In Fourier space, it turns out that all self-energy insertions resum in such a 
way that the equilibrium contribution attains a finite width, while the out-of-equilibrium 



terms cancel 33| . We briefly review this in Section [2j It is sometimes argued that the 



occurrence of pinch singularities signals a breakdown of perturbation theory for out-of- 



equilibrium systems, even when the coupling constant is arbitrarily small (3J, |35[. We 
emphasise however, that within the Schwinger-Dyson approach, no such problem occurs. 

The zero-width behaviour also appears to be related to a feature of the Schwinger- 
Dyson equations on the CTP that is discussed in Refs. 0, [32[. When interpreting the 
kinetic equations f llSbp . that appear as a subset of the Schwinger-Dyson equations, as 
Langevin equations, they obviously feature a Stokes term, that induces the approach of 
the macroscopic observables to equilibrium, but they lack a Langevin term, that induces 
fluctuations and might lead to a finite width for the non-equilibrium distribution. In 
Ref. 32|, the kinetic equations are reformulated in a way such that they take a form 
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similar to Langevin equations. Since the resulting equations are however equivalent 
to the original kinetic equations in the stationary bath approximation, the feature of 
zero width distributions still persists within the solutions. The truncation of n-point 
functions is proposed as the reason for the absence of a Langevin term in Ref. 0]. 
Clearly, when the effects of backreaction are neglected, the functional dependence of the 
self-energies on the two-point functions is dropped in order to achieve a linearisation of 
the Schwinger- Dyson equations. However, while stochastic sources for the fluctuations 
of the propagators ~ ~ addressed, the zero width contributions to A^'^ 

are not discussed in Ref. j7|. 

The occurrence of zero-width terms appears counter-intuitive. Based on the expecta- 
tion that one should not be able to tell whether a quasi-particle is part of the equilibrium 
or non-equilibrium distribution by measuring its width, rather than the solution f lT8|l . 
one might use the ansatz f l22|) . where all quasi-particles have a finite width. Naively, 
we might expect important consequences for calculations, that rely on the heuristic 
ansatz fl2^ . where both the equilibrium and the non-equilibrium contributions to the 



quasi-particle distributions have a finite width, as e.g. in Refs. |16H21l. I36l . I37j . where 
off-shell contributions to the scattering rates are of leading importance. 

A purpose of the present work is to resolve this seeming contradiction. In Section |2l 
we review the summation of pinch singularities in the CTP formalism. The propagator 
with zero-width out-of-equilibrium contributions in the limit of negligible backreaction 
is derived in Section 121 An explicit numerical check of these solutions for a scalar theory 
with quartic interactions and a single excited mode is performed in Section |H The 
Figures in that Section illustrate the behaviour of the non-equilibrium and equilibrium 
components of the statistical propagator (the sum of the two Wightman functions) in 
two-time representation. These discussions confirm the validity of the solution flTH]) for 
the Wightman function with a zero-width out of equilibrium component in situations 
with negligible backreaction. The agreement between numerical and analytical results is 
of great value, because we can exclude that the zero-width behaviour is due to analytical 
approximations that are not understood. In Section [5l we show that the summation 
of all gradient terms, that act on the on-shell 5-function leads to an effective finite- 
width behaviour of the out-of-equilibrium distributions in the collision term. In order 
to "observe" the finite width, in Section [6l we propose to use interactions that are 
kinematically forbidden for zero-width particles. The main conclusion of Sections |5] 
and [6] is that the use of the solution f|T8l) with the zero-width components and the 
summation of the gradient terms is equivalent to the use of the heuristic ansatz fl22l) 
without summation of the gradient terms. This way, the solution f[TK]l is reconciled with 
heuristic approaches, that assume a finite width for the out-of-equilibrium distributions. 
We may also interpret this result as a generalisation of the fluctuation-dissipation relation 
to out-of-equilibrium systems in situations, when backreaction is negligible. Further 
discussions and conclusions are offered in Section [71 
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2 Pinch Singularities 



The analysis of this Section is originally reported in Ref. [33[. Here, we recapitulate it in 
order to put it more directly within the present context. We consider a real scalar field 
(yj of mass m^p within the CTP formalism. In this Section, we evaluate the propagator, 
which is a two-point function, in Fourier (momentum) space. As an implicit consequence, 
we do not allow for a dependence of the propagator on the average coordinate, i. e. we 
cannot expect it to describe the approach of the field ip toward equilibrium. The tree-level 
CTP propagator for the scalar field can then be expressed in matrix form as 



iA(°)^(fc) 




Uk)U{k)\_.^^,^^ f u{k) uik) 

Uk) Uk) ; ^ ^ ^ I l + /^(A;) l + ^(A;) 



where 



:::: 

and /(p(k) is the distribution function in momentum space, and 

are the tree-level retarded and advanced propagators, with e being infinitesimal. For the 
definition of the various two-point functions on the CTP and their relations amongst one 
another, cf. Appendix |X1 The tree-level Wightman functions have zero width: 

\^f<{k) = iAW+-(A;) = 27i5{k^ - ml)sign{k')U{k) (4a) 

= 27r5(A;2 - mj) [^(A;o)^(k) + ^9(-fco)(l + ^(-k))] , 
iAf>(A;) = iAW-+(A;) = 2n6{e - ml)sign{k')[l + U{k)] (4b) 

= 27T6{e - ml) [^9(A;o)(l + /^(k)) + ^(-fco)^(-k)] . 
We note that the zero-width spectral function 

2Aj')-^(fc) = 27r5(fc2 - mysign(/t°) (5) 

occurs as a factor within these expressions. 

We aim to obtain Wightman functions of finite width by summing over all insertions 
of self energies ill^, which are understood to be CTP matrices. For this purpose, we 
define the matrix 



(6) 



to take account account of factors ±1 that occur in the CTP Feynman rules. The self- 
energies in^ on the CTP are defined in analogy with Eq. ([1]) . For insertions in the 
tree propagator, we obtain 

iA(°) (sn^siAj))^ (7) 
= i A(°)^ (m^Aj')^) ""('^yA- i A?)^ (in^iAj')^) ^ ^ ^' 



iA(°)^iA(°)^ [Um> - (1 + u)m<] ^(iA(°)^inJ)^-(iA(°)^in 



Af-l 
j=0 



1 1 



where we have suppressed the arguments k. The last term contains the so-called pinch 
singularities, since products of poles with the same real parts but opposite imaginary 
parts occur there. From this expression, we obtain the resummed, full propagator 



iA, = 5^iA(°)(Sn,SA(°))^ (8) 



N=0 

^ i f ^ + U U\ i f U U 

-i^i^ [«-...«,(;; 

where the last term results from the resummation of the pinch singularities. When 
we make use of = i (n> - n<) = i (n^ - H^) (c/. Appendix lA]) and define the 
stationary distribution function 

f-Uk) = ^^^^^ (9) 

the resummed propagator is given by 

iA (k) ' ( ' + ^^^^^ f-^^^ ] m 



k'-ml-U^ yi + f^^(k)l + f^\k) 
In particular, the Wightman functions are now 

- m2 - my + 



iA<(fc) = iA+-(A:) = — . Z^.^ . = 2A^(A;)/-(fc) , (11a) 



211-^ 
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which have the desired finite width for 7^ 0. But evidently, there is the unwanted 
feature that all dependence on the distribution function /^(/c) has vanished in a spec- 
tacular way. Instead, the propagator depends on f^{k), which is determined by the 
self-energy Il^{k). This is a consequence of the fact that the calculation has been per- 
formed in Fourier space, with no provision for a time dependence that would describe 
the relaxation of the system toward equilibrium. The only possible solution is the one 
that is in equilibrium with the "background" described by the self-energy H,^. 

The Schwinger- Dyson equations in Wigner space in principle account for the time 
dependence. In Section we show that in order to capture the finite- width behaviour for 
the time-dependent out-of-equilibrium contributions, one must include terms of all orders 
in the gradient expansion as well. The finite width of the out-of-equilibrium distribution 
appears not as a feature of the propagator itself, but it is rather encompassed in the 
collision term expanded to all orders in the gradient expansion. 

3 Solutions to the KadanofF-Baym Equations with- 
out Backreaction 



In this Section, we solve the Schwinger Dyson Equations on the CTP under the assump- 
tion that all self-energies Ilip{x,y) are time-translation invariant. (In Wigner space, this 
implies that they are time-independent). This amounts to an approximation where the 
backreaction of the field if on the fields that are relevant for the self-energies 11^ is ne- 
glected. We derive the solutions in Wigner space and thereby recover the result that has 
been found earlier in two-time representation within Ref. 

For a real scalar field ip of constant mass m^, the Kadanoff-Baym equations in Wigner 
space are 



- -dl + ik'^dt - ml 



e-'^{nJ}{A<'>}-e-^{n<'>}{An (12) 



= 2«-^M{n>}{A<}-{n<}{A>}) 

and the equations for the retarded and advanced propagators 
1 



F + ik^'dt - -df 
4 



m,, 



^R,A _ eio|nJ}{Af ^} ± e-^{in^}A«'^ = 1 , 



where the diamond operator is defined as 

1 



o{A}{B} = - (d^A) {dkB) - - (dkA) (d^B) 



(13) 



(14) 



and where t = x^. Provided the scale of macroscopic variations ~ is smaller than the 
typical energy-scale ~ {dk)~^, we may expect that we can truncate at a finite order in the 
diamond operator. This approximation scheme is known as the gradient expansion. We 
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refer to the right-hand side of Eq. f|T2|) as the colhsion term. In Section [5l we note that the 
identification of ~ {dk)~^ with the energy scale does not hold, when the derivative acts 
on an on-shell (5-function, but also that these terms can be straightforwardly summed 
and interpreted. 

It is useful to split the Kadanoff-Baym equations into constraint and kinetic equations 

as 



A<'> - cos(o){nf }{A<'>} - cos(o){n<'>}{Aj} (15a) 

= -^sin(o) ({n>}{A<} - {n<}{A>}) , 

A;°9,iA<'> + sin(o){nJ}{A<'>} + sin(o){n<'>}{Aj} (15b) 

= -^cos(o) ({in>}{iA<} - {in<}{iA>}) . 

The assumption of a constant 11^ amounts to neglecting its functional dependence on 
A^, and Eqs. fll2|13p become linear, such that they can be solved using straightforward 
methods. Of course, this linearisation corresponds to the neglection of backreaction. 
In particular, the result for the stationary propagator from the resummation of pinch 
singularities is reproduced easily: Eq. ( IT3|) yields 

= fc^-m^-n^±in^ ' 

A^ = -rA^ + A^) = ^ ^ (16h) 

which can be substituted in Eq. fll5ap such that 

(P_^2 _nH)2 + n^2 .^n>-n<) -mi-n^y + uf^^ ' 

(17) 

where we define f^^ = and f^^ = 1 + f^. The resummation of pinch singularities is 
implicit within the Schwinger-Dyson equations on the Closed Time Path, such that their 
solution is perhaps a simpler and more direct procedure for obtaining approximations to 
the full propagators. 

Under these assumptions, the Wigner space solution to Eqs. (fT^ and ( IT5]) at zeroth 
order in the operator o is 

iA<'>{k,t) =iA;*<'>(A;) + 6U{k',k)sign{k')27T6{e - - n^)e-(n^/'=")* (18) 
=iA;,*<'>(A;) + 6U{k)27r6{e - - n^)e-r.(k)* , 
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where (5/^(A;°,k) is an integration constant. Note that the solution is always decaying 
for both particle (/c° > 0) and antiparticles {k^ < 0), because Il';^{k^,k) is odd in /c°. 
The obvious interpretation of 6f^(k) is the deviation of the particle distribution function 
from its equilibrium form. 

For the second equality, we have used the Breit-Wigner approximation. This corre- 
sponds to the assumption that it is sufficient to know the value of li-^ik) at the quasi- 
particle pole, where 



= ±u^ik) = ±y^k2 + ml + n^K(k), k) . 

The approximation is valid when only contributions close to the pole are important and 
n^(/c) is a smooth function in k^. In other words, we may approximate 

n^(At, k) = la;J(k)r^(k) sign(At) , (19) 
where uj^{k) = -^/k^ + . Then, 

n^(^°'k) = ^r,(k) (20) 

is approximated by a function which is smooth in k'^ close to the quasi-particle pole. Here, 
r<^(k) is the relaxation rate. In the present context, this approximation is obviously valid, 
because the ^-function forces the evaluation of U.'^{k) at the quasi-particle poles. 

The fact that 11^, which appears within Eq. (fTSl) as the finite width for the stationary 
(equilibrium) distribution f^, and therefore describes fiuctuations, is the same as in the 
retarded propagator fll6ap . which describes the dissipation of a linear perturbation to- 
ward equilibrium, is a consequence of the celebrated fiuctuation-dissipation relation for 
equilibrium systems. In this light, the absence of fluctuations in the terms that multiply 
6f^ in Eq. f|T8|) may appear as a surprise and suggest that the fluctuation-dissipation 
relation cannot at all be generalised to out-of-equilibrium situations. However, in Sec- 
tional we show that even though the non-equilibrium contribution to iA|^'^ in Eq. flTBl) 
does not exhibit a width, fluctuations arise effectively from a summation of the gradients 
within the collision term. 

The inverse Wigner transform of Eq. (fTSi) then leads to the result of Refs. 11, [l3|l^ 



1 ry,(k)sign(At) 



iA<'>(t;At)^-e^ ^^T— At 



(21) 



^ The solutions in Refs. ll|, ll3| also encompass possible oscillatory terms in the average coordinate 
t with the frequency 2u!. These can be identified with the coherence solutions that are reported in 
Refs. [H-Iill- 
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where we have expanded to hnear order in 11;^ within the exponentials and the prefac- 
tors. Apparently, the distribution f|T8|) consists of a finite-width equilibrium term and a 
zero-width non-equilibrium contribution. It naively appears that we could tell whether a 
particle belongs to the equilibrium part or the non-equilibrium part 5f by observing 
its width. In the coordinate representation (ETj) . this is refiected by the fact that the 
correlations proportional to the out-of equilibrium distribution Sf exhibit no damping 
in the relative time coordinate At. A corresponding behaviour is exhibited by the prop- 
agators for out-of-equilibrium fermions reported in Refs. 
the same assumptions. 



121, Q, that are derived under 



The solution f ll8p should be contrasted with the "Kadanoff-Baym ansatz". 



iA<(A;) 



2U^ 



2 fv(k) 



2n^ 



:[l + U{k)]=2A^ik)[l + U{k)]. 



(22a) 



(22b) 



This relies on the intuitive expectation that it should not be possible to distinguish be- 
tween an equilibrium and a non-equilibrium particle by measuring its width. In Section |6l 
we describe how such a measurement may be performed. This ansatz is proposed more 
or less explicitly within the literature on Kadanoff-Baym equations, see e.g. Eq. (4.44) 
in Ref. 32|, Eq. (11.13) in Ref. 15| or Refs. jsl, [s^j- In principle, it also underlies cal- 
culations in more heuristic app roaches to kinetic theory, where finite-width effects are 



important, such as in Refs. (I6l-l2l[|. In the subsequent Section HJ we numerically confirm 
the correctness of the solution ( IT8|) . Nonetheless, as we show in Section [5l it turns out 
that the ansatz (122|) effectively resums an infinite number of gradients within the collision 
term, such that its use in kinetic theory can be justified. 



4 Out-of-Equilibrium Propagator from a Numerical 
Simulation 

In this Section, we present a numerical check that indeed, when backreaction effects may 
be neglected, Eq. ( IT8|) is a correct solution to the Kadanoff-Baym equations, while the 
naive finite- width ansatz ( 122|) is not. 

The assumption that is crucial for deriving the analytical solution ( !T8|2T]) to the 
Kadanoff-Baym equations is the time-independence of the self-energy in Wigner space, 
n<^(t, /c/i) = Il^{kfj). This situation is realized for example when the non-equilibrium 
field (p interacts with a thermal bath. In this case, a necessary condition for the time- 
independence of the self-energies is that the thermal bath can be considered as an in- 
finitely large reservoir, that is not infiuenced itself by the interaction with the non- 
equilibrium field ip. 

In this context the question arises how sensitive the analytical solution (]18|2ip . and 
in particular the peculiar behaviour of the non-equilibrium contribution, is to the under- 
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Figure 1: Statistical propagator A;^ = (A^ + Aj^)/2 (red) and spectral function A;^ = 
^(A^ — ^'^) (blue) obtained from a numerical solution of the Kadanoff-Baym equations. 
Shown is the dependence on the relative time At for fixed central time t. 



lying thermal-bath assumption. In order to investigate this question, we compare the 
analytical expression with numerical solutions to the complete Kadanoff-Baym equa- 
tions, for which the self-energies are computed self-consistently. This ensures that the 
full back- reaction is taken into account. Concretely, we consider a real scalar field with 
Xip^ interaction, and solve the Kadanoff-Baym equations in the two-time representation. 



+ m, 



X 



-A<'>i 
2 



x,y) 



to 



d'zU<'>ix,z)A^{z,y) 



2 r d'zU^{x,z)A<'>{z,y) 

J to 



where to is the initial time, A;^ = ^(A^ — A;^) is the spectral function, and 



d^z 



to 



to 



dz^ / d-^z 



For the self-energies, we use the so-called setting-sun approximation. 



n 



x,y) = -yA<'>(x,?/)' 



(23) 



(24) 
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relative time At [1/mth] 



Figure 2: As figured! but for the excited momentum mode |k| = S.Omth. The statistical 
propagator = (A^ + Aj^)/2 can be described by the sum of an exponentially damped 
equilibrium contribution and an undamped non-equilibrium contribution. 



which describes scattering processes, and is accurate for small enough coupling and 



occupation numbers (see [42|, |43| for details). We assume that the system is spatially ho- 
mogeneous and isotropic. In this case it is convenient to transform to spatial momentum 
space, 

(25) 



A<'>(t,At,|k|) = j d'xe-''^<''-^^A<'>{x,y) 



where, in addition, we have introduced the central time t 



{x^ + y^)/2 and the relative 
time At = x° — for convenience. The Wigner space functions would result by an 
additional Fourier transformation with respect to At. 

Since we are interested in situations where a non-equilibrium degree of freedom inter- 
acts with many other degrees of freedom that are close to thermal equilibrium, we initial- 
ize the system at time t = At = with thermal correlators at temperature T = l.brrith, 

where rrith is the thermal mass. As described in i ^ 

Gaussian four-point correlation into account. In addition, 
to be |k| = S.Orrith, is excited above the thermal occupation with an effective particle 
number An = 7.5 (see 42|). For the numerical solution we use a lattice of size 1400^ x 32'^ 
in (a;°, |/°, k)-space, a time-step at = 0.025/mth, and momentum cut-off A = n/a^ with 
Qx = 0.5/mth- Furthermore, we choose Ab/4! = 0.75 for the renormalized coupling and 
determine the counterterms according to |45[. The numerical solutions conserve the to- 
tal energy with an accuracy of better than 0.05% (note that this property would be lost 



we also take the thermal non- 
one momentum mode, taken 
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when neglecting the back- reaction) . 

The plots show the dependence of the two-point functions on the relative time At for 
the excited momentum mode |k| = S.Omth, as well as a mode with slightly smaller 
momentum |k| = 2.96mf/j. All correlators are evaluated at the fixed central time 
t = 17. 5 /rrith- For the non-excited mode shown in Figure [H we find that both the statis- 
tical propagator and the spectral function are damped exponentially with respect to the 
relative time, as expected. Furthermore, the damping rate is identical for both correlation 
functions. This property is expected from the fluctuation-dissipation relation which is 
valid in thermal equilibrium. In addition, the relative magnitude of the statistical propa- 
gator and the spectral function is close to the expected value l + 2fBEi^/^^ + T^th) — -'-•3' 
where Jbe is the Bose-Einstein distribution function. 

In Figure [2l the corresponding propagators are shown for the excited momentum 
mode |k| = S.Orrith- The spectral function is also damped exponentially, whereas the 
statistical propagator shows a different behaviour. However, the latter can be described 
to a good accuracy by the sum of the exponentially damped equilibrium contribution and 
an undamped non-equilibrium contribution. This is precisely the behaviour predicted 
by the analytical solution fll8|2ip . In addition, we have checked that when varying the 
central time t, the non-equilibrium contribution to the statistical propagator decays ex- 
ponentially, while the equilibrium contribution remains constant. In addition, the decay 
rate of the non-equilibrium contribution obtained in this way coincides with the damping 
rate of the equilibrium propagators. Again, this result is in accordance with the analyti- 
cal solution f fT8|2T]) . We conclude that the analytical solution obtained under the thermal 
bath assumption constitutes a reasonable approximation to the full numerical solutions, 
and is in particular consistent with a non-equilibrium correlator that is undamped with 
respect to the relative time. 

5 Effectively Finite Width from the Gradient Ex- 
pansion 

We have now established analytically and numerically, that when backreaction may 
be neglected, the Wightman functions ( ITS]) that feature zero-width out-of-equilibrium 
contributions are indeed correct solutions to the Schwinger-Dyson equations on the CTP. 
In this Section, we investigate, whether the finite- width ansatz fl22|) for iA^''*, that 
we might intuitively expect and that it is often assumed in the literature, may yet be 
appropriate, and if yes, in what sense. 

For this purpose, we substitute Eq. (fT8|) into the right hand side of Eq. ( |T2ll . The 
substitution results in 

J2 - ( -^r^k) [(4o)"n^] SU{k\ k)sign(A;°)2vr5(A;2 - - nJ)e-^-«* . (26) 

n=0 V / 

We now use the fact that we evaluate these terms in the distributional sense, after inte- 
gration over such that we can swap the /c°-derivatives through integration by parts. 
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Therefore, within the (5-functions, we replace — k'^+^T^, because Eq. fl26|) corresponds 
to a Taylor series. One may worry about the fact that we have deliberately replaced 
expressions of n^//c° by r^(k) and therefore truncated additional dependence. How- 
ever, we may estimate derivatives with respect to acting on these factors as yielding 
factors of T~^, such that the error incurred within each order n of the gradient expansion 
can be estimated as (r(^(k)/T)". The same argument applies to derivatives acting on 
5f^{kP, k). Derivatives acting on sign(/c°) are suppressed because — — 11^) = 
for = 0, at least provided + 11^ ^ 0. The suppression of gradients does obviously 
not apply to derivatives acting on the (5-function. The consequences of this is what we 
explicitly calculate in the following. 

From Eq. fl2Tl) . we see that the (5-function originates from the Fourier transform of 
cos(a;^(k)At) with respect to At. The real part of Eq. fl26l) enters the kinetic equa- 
tions fll5bp . while the imaginary part enters the constraint equations fllSap . To this end, 
k^ has been assumed to be real, while now, we aim to continue the expressions to com- 
plex k^. This implies, that we need to identify the various occurrences with k^ or k^* . 
The correct prescriptions are determined by the causal properties of the propagators, 
i. e. in two-time representation, whether the time argument is evaluated above or below 
the real axis, or in Wigner space, where the poles are located with respect to the real 
axis of the /c°-plane, cf. Appendix [B] for a quick reminder on this matter. In two-time 
representation, the correct prescription for the imaginary part of the time argument is 
given by 

oo 

2A4(A;)=^i^^ j rfAte(^'="-^^-(''>^^'^(^*»^*cosKAt) (27) 



sign(A;°) 
" 2u;^(k) 



E 

± 



A:0±u;^(k) + ^r^(k) A:0±a;^(k)-^r^(k), 



sign(fc°) r^(k)(fc°^ + u;^(k)) 

" a;^(k) _ ^2(k))2 + po + ^^(k))2 + _ ^2(k))2] ir2 (k) + ±r^(k) 



(F2-a;2(k))2 + p2r2(k) 

This result is purely real, such that it enters into the kinetic equations, but not into the 
constraint equations. Note that for these equations to be valid the actual integration 
interval must be large compared to F"^. Regarding the replacement k^"^ — )■ c<;^(k) in the 
numerator, one may worry that one drops a term that gives rise to relevant contributions 
to the integrations over k^. Note however, that within the collision term, this finite width 
representation is always multiplied by distribution functions. We therefore assume, that 
all distribution functions decay exponentially for large values of as it is the case for 
e.g. the equilibrium distribution functions. 
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We also observe that 



(k) 



1 



(28a) 



m 



2 + iefco ' 




1 



) 



* 



A;2 — + iek^ 



(28b) 



may be continued to complex k^ with lm[A;'^] > while maintaining the causal properties, 
because these expressions are analytic in the upper complex half-plane. In terms of these 
propagators, we may express 



For real k^ and infinitesimal e, this corresponds to the representation of 27T6{k'^ — — 
n;^(^)) = 2A'^°)-^ as in Eq. fl26l) . Through the Taylor expansion, this expression can be 
directly continued to complex arguments as — )■ + jT^{k) 



in agreement with Eq. f l27|) . 

It should be noted that due to the singular behaviour of the on-shell (5-functions, it is 
not consistent to truncate the gradient expansion of the collision term at any finite order 
in the diamond operator. Here, we have shown that the gradients can easily be resummed 
to yield a benign contribution with an intuitively clear meaning: an effective finite-width 
distribution function. We therefore expect that the systematics of the gradient expansion 
still applies to kinetic theory, once singular terms arising from on-shell 5-functions are 
summed, following the methods introduced in the present Section. 

In conclusion, the higher order contributions from the gradient expansion that arise 
within the collision term precisely recover the naive finite- width ansatz fl22|) ! In this 
sense, the fluctuation-dissipation relation generalises to an out-of-equilibrium system of 
the scalar field (p, for which backreaction can be neglected. In turn, when we substitute 
the finite width ansatz (l22l) into the collision term, we must not expand to higher orders 
in gradients, as this is already accounted for by the ansatz itself. 



6 Application to Processes with Kinematic Thres- 



Provided the self-energy Il^{k) is not changing rapidly in k^ close to the pole, there 
is no leading order difference in whether we integrate over the (5-function or the finite 
width representation of the propagator. This behaviour is however precisely what we 
assume within the Breit-Wigner approximation (fT9l) . Moreover, when we integrate the 



2A-^ 



P - m2 - m (k) + iek^ (F)* -ml- U^{k) - iek^ ' 



(29) 




(30) 



holds 
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kinetic equation (I15b|) for from — oo to oo, both sides are trivially zero, because of the 



symmetry properties of the real field. Hence, in order to "measure" the effectively finite 
width behaviour found in the previous Section, we yet need to specify an interaction 
that can probe the off-shell contributions of iA|^'^(/c) for ^ a;^(k). 

We therefore consider a reaction that is mediated by a coupling h and that is kine- 
matically forbidden at tree-level. For such a process, accounting for the effective finite 
width of the propagators is indeed of leading order importance. Second, in order to ex- 
tract non-trivial information about the kinetic evolution, instead of plainly integrating 
over k'^, one needs to include integrals over higher moments (additional powers of k^ 
within the integrand). For example, the distribution function can be obtained from the 
integral 

/.(k) = / ^|fc°|iA<(A:)-l. (31) 

For different observables, analogous prefactors for iA^ can be found that are typically 
smooth functions of k^ in the vicinity of the quasi-particle pole, such that derivatives 
may be consistently neglected at leading order in the gradient expansion. 

For the purpose of probing the width of the field ip, we introduce two additional 
scalar fields xi,2- The equation that determines the evolution of /^^ follows from the 
kinetic equation (]15bp (with the obvious replacement f ^ Xi) 

dtfA^)=j ^l^°l^*i^x;" (32) 

— sign(A:°)-cos(o) ({in> }{iA< } - {in< }{iA> }) . 

Here, we have truncated terms that are of higher order in gradients on the left-hand side. 

In the presence of kinematic thresholds, including the gradient terms that lead to the 
effective finite width can be of leading importance. In order to illustrate this, we consider 
an interaction described by the Lagrangian term hipxiX2, where /i is a real dimensionful 
coupling constant. For the field Xi? this induces the Wightman-type self-energies 

ill<i^{k) = h' I ^iA<'>(A:-p)iA<f (p). (33) 

Suppose now, that the field ip deviates from equilibrium, as described by the distribution 
function Sf^{p), while the fields Xi,2 are in equilibrium at the temperature T. The fields 
have the masses m^p and m^^ ^ . Moreover, assume that the production or decay rates of 
Xi,2 due to the interaction with ip are much smaller than the decay rate F^. Then, we 
can use for iA^^ ^ the on-shell propagators. Moreover, the leading time dependence is in 
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(5/<^(p), and we may evaluate the action of the diamond operators as 
1 



-h'Re 



oo 
I- n=0 

i 



ni 



d^p rdk^ r dp^ 

(27r)3 J ^ J ^ 
iA>ik-p))d^,oiA<^{k)signik') 



(34) 



We could now directly apply the Taylor expansion to the propagator A^-^. However, 
since it is the field (p which is the origin of the effective finite width behaviour, we rather 
perform an integration by parts, in order to attach the derivatives dk^ to iA^p. We again 
ignore derivatives acting on T^, since these would lead to terms that are suppressed by 
higher orders in the gradient expansion. Note also that we can ignore the derivatives 
acting on the sign function, since the on-shell propagator iA^^(fc) vanishes for = 0. 

Substituting the explicit forms of the propagator and performing the Taylor expansion 
of iA^{k — p) toward complex k^ , as it is explained in the previous Section, we obtain 



X 



d^p 
2r^(k- 



dk^ 
27 



27 
-P) 



sign(p°)27r(5(|}^ - m'^_j2nS{k 



m 



(35) 



[{k - pY - mlf + {k^^ - p^fVl 



(^/■^(k- p)e 



-r^(k-p)t 



.kO/T 



Now assume for definiteness, that m^-^ + m^^ > but m^^j + > m^^-^. This is 
for example the case when all masses are of equal or similar size. A 1 -H- 2 process 
between ip and Xi,2 then is not possible provided all particles are on-shell. In the present 
case, the field (p obtains an effective finite width, and dtf^iO^) is non-zero, even for 
the given mass relations. The above integral can in general be evaluated numerically 
in a straightforward way. Let us however additionally assume that m^j2,m^ ^ T and 



m 



XI 



+ m 



X2 



Xl,2,f 



(k) 



-C m^j2,m^. In that situation, all particles are non-relativistic such 



m 



Xl,2,^ 



and r^(k) r^(0). Moreover, the dominant reaction channel is 



that uj- 

X1+X2 ^1 where ip is off-shell. The meaning of ip being off shell is that it corresponds to 
a virtual state, that eventually rescatters with addtional particles, which is encompassed 
within the self-energy 11^. In this situation, the angular integration can be performed 
easily. Moreover, we note that in the integrand sign(p'^) 7^ sign(/c''). With the above 
approximations, we therefore obtain 



1 l2 1 [ P^dp 1 



4r^(0)m, 



2 2a;^,(k)7 (27r)2 2u;;,,(p) r^2 , 

L ^1 

1 

X y"rfcos^9 5/^(k-p)e-^^(°)*. 



(36) 
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where i? is the angle between k and p and where we have assumed that 6f^(k) is isotropic. 
Clearly, this can be interpreted as a scattering process, between ip and additional par- 
ticles, that are described by F^, via a virtual off-shell ip into Xi X2- In contrast, 
the result would be strictly zero if we would not have taken account of the resummation 
leading to the effective finite width propagator as discussed in the previous Section. We 
emphasise that for the present calculation, it is sufficient to take the Breit-Wigner ap- 
proximation for n^(/c), while the dependence of 11^'-^ (k) on is of leading importance. 
For related discussions of off-shell effects, cf. Refs. [13|, |37 . 



7 Conclusions 

For an out-of equilibrium scalar field, that interacts with a time-independent, spatially 
homogeneous bath and that only experiences negligible backreaction, we have reviewed 
and derived various aspects of the solutions to the Schwinger- Dyson equations on the 
CTP. Our main focus is the solution for the Wightman function (ITSll and its curious 
feature, that it decomposes into a finite-width equilibrium and a zero-width out-of- 
equilibrium contribution. The corresponding Wightman function has been derived before 



in position space in Ref. but no direct interpretation of the different width of the two 
components is provided there. The fact that the equilibrium component exhibits a finite 
width is confirmed by a direct summation of the geometric series of insertions of the 



self-energies into the propagators. This calculation has first been performed in Ref [33 
and in Section [2] we have reviewed it in order to put it within the present context. The 
desirable cancellation of pinch singularities is accompanied by a spectacular cancellation 
of the out-of-equilibrium contributions. In order to verify that the vanishing width of 
the out-of equilibrium component is not an artefact of the analytical approximations, 
in Section HI we have presented a numerical study that clearly confirms this feature. 
We emphasise that when using the Schwinger-Dyson approach, we do not encounter a 
problem with pinch singularities. Besides, we do not find a breakdown of perturbation 
theory applied to out-of-equilibrium systems, that is sometimes associated with pinch 



singularities [34, l35 



The zero width behaviour of the out-of-equilibrium Wightman function, that we have 
carefully confirmed calls for interpretation. Intuitively, we expect that it should not be 
distinguishable, whether a quasi-particle is part of the equilibrium or out-of-equilibrium 
distribution. In Section |5l we present the observation, that an effective finite width 
for the out-of-equilibrium contribution results from a summation of all gradients within 
the collision term. Such a summation to all orders in gradients becomes necessary, 
because derivatives acting on the on-shell (5-function are not suppressed and therefore 
do not comply with the usual suppression of higher orders derivatives in the gradient 
expansion. Section [6] provides an example, how the effective finite width is observable 
through interactions, that are kinematically forbidden for zero-width particles. Since the 
effective finite width takes the same form as the width of the equilibrium propagator, 
our findings may be interpreted as a generalisation of the fluctuation-dissipation relation 
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to out-of-equilibrium systems with negligible backreaction. 

The results of the present work provide a justification for the use of a heuristic finite- 



width ansatz, that is often intuitively applied within kinetic theory, e.g. in Refs. |16H21 



36|,l37[. If one aims to derive kinetic theory from first principles of Quantum Field Theory, 
using the CTP formalism, the present work fills in a missing step in order to describe 
the finite width of out-of-equilibrium particles. Note that the propagators for fermions 



reported in Refs. [12|, also feature a zero-width non-equilibrium contribution. The 



resummation technique leading to an effective finite width can in principle be generalised 
to fermionic systems. 

In the field of leptogenesis, we see future applications in studying the role of the finite 



width of the singlet neutrino, a discussion which has been initiated in Refs. |12l . 114 . 
Moreover, the correct description of the finite width of out-of-equilibrium neutrinos is 
potentially crucial for resonant leptogenesis 46|-|50| in a regime, where the poles of the 
nearly mass-degenerate neutrinos overlap within their width. 

In conclusion, we have shown that the Schwinger- Dyson equations on the CTP in 
the linear approximation, that allows for analytical solutions, are well suited to describe 
finite-width effects for out-of-equilibrium systems. Due to the zero-width behaviour of 
the Wightman functions, this is not directly obvious, and to arrive at this result, a non- 
trivial summation of gradients of all orders acting on on-shell ^-functions is necessary. 
Once this summation is performed, it appears that the usual gradient expansion can be 
systematically extended order by order. In the future, we may therefore expect further 
progress in the application of the CTP formalism to the analytical description of out-of- 
equihbrium systems. 
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Note added: During the preparation of this Article, Ref. [51j appeared, which contains 
a similar argument about an effective finite width, that does not arise from the relaxation 
toward equilibrium, but from a time-dependent mass term. 



A Two Point Functions on the CTP and Wigner 
Space 

Within the present work, we follow the notations and conventions as given in Ref. joj. 
For a quick reference, we quote here the relations between the two-point functions, that 
are used within the present work, where G may either stand for a CTP Green function 
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A or a self-energy 11: 

= - G'> = G'< - (advanced) , (37a) 
G^ = G^ -G< = G> - (retarded) , (37b) 

G^ = i(G^ + G^) (Hermitian) , (37c) 

G^ = ^(G^-G^) = ^(G>-G<) (anti-Hermitian, spectral) . (37d) 

The identification of these two-point functions with those bearing CTP indices is 

= G++ , G< = G+- , G> = G-+ , G^ = G— . (38) 

The functions G^'^ are known as the Wightman functions, G'^ is time ordered, G^ 
anti-time ordered. 

The Wigner transform of a two-point function is defined as 

G{k, x) = j dW''G{x + r/2, X - r/2) , (39) 

where we refer to r as the relative and to x as the average coordinate. In the situation of 
spatial homogeneity, there is no dependence on x (spatial translation invariance), such 
that G{k,x) = G{k,t). 

B Prescription for the Imaginary Part of the Time 
Argument 

The time-ordered and anti-time ordered Green functions of a scalar field (p are 

iA^'^Ax) = / T^J^e--^^^^, (40) 

J (27r)4 pO -u^{p)±ie 



0\ loip-Ax 



J (27r)32a;^(p) J 

Thus, the prescription for the imaginary part always appears in such a way that the inte- 
grand decays exponentially for large lAx*^]. Also recall that we may write iA'^''^(Aa;) = 
'?9(Aa;)iA^'^(Ax) + ??(— A,x)iA"^'^(Ax). This explains the explicit appearance of sign(At) 
within Eq. ([2] 
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